LARGE TIME BEHAVIOR OF THE HEAT KERNEL OF TWO-DIMENSIONAL 
MAGNETIC SCHRODINGER OPERATORS 

HYNEK KOVARIK 

Abstract. We study the heat semigroup generated by two-dimensional Schrodinger operators with 
compactly supported magnetic field. We show that if the field is radial, then the large time behavior 
of the associated heat kernel is determined by its total flux. We also establish some on-diagonal heat 
kernel estimates and discuss their applications for solutions to the heat equation. An exact formula 
for the heat kernel, and for its large time asymptotic, is derived in the case of the Aharonov-Bohm 
Cn , magnetic field. 
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1. Introduction 



The Hamiltonian of a charged quantum particle in R d interacting with a magnetic field B is given 
formally by the differential operator 

H B = {iV + A) 2 (1.1) 

in L 2 (M. d ), where A is the vector potential of the magnetic field; B = rotA (for d = 2,3). The 
object of our interest in the present paper is the integral kernel e~ tHB (x,y) of the heat semigroup 
generated by Hg. In particular, we are interested in the dependence of e~ tHn (x,y) on t. A well 
known semiclassical result, [Er97, Ma] says that under certain conditions on B we have 

lim ti e- tHB (x, x) = (47r)- d/2 . (1.2) 

In other words, the leading term of e (Hb (x, x) in the short time limit is not affected by the magnetic 
field. However, the situation changes in the large time limit, where the diagonal element of the heat 
kernel decays exponentially fast provided the size of the magnetic field is bounded from below by a 
positive constant, [Er94, Mai . More precisely, the estimate 

lim - log||e- tHlB |Ui_>£oc < ~C L L min ||B(a;)|| (1.3) 



holds true with Cl = 1 + o(l) as I - ► oo, see |Er94) . From the Mehler formula for the heat kernel 
of the two-dimensional Schrodinger operator with a constant magnetic field, see [Sil| Sect. II. 4-6], 
it follows that the factor min^g^d ||B(x)|| in (|1.3|) cannot be improved. Later, a uniform pointwise 
bound on the two-dimensional magnetic heat kernel in the form 



■"';;il. .. r ^L 

47rsinh(^) ' '' x6M 2 



e- tHB \\ L i^ L ~ < - , /' -T- , B = min \B(x) |, t > 0, d = 2 (1.4) 



was obtained in jLT) under the assumption that Bq > 0. This bound is the best possible since there 
is equality for B = Bq. The latter follows again by the Mehler formula. 

In this paper we focus on the the case d — 2 and address the following question: what is the large 
time behavior of e~ tHB (x, y) when B{x) is of compact support? Note that for a compactly supported 
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magnetic field we have -Bo = min^g^ |-B(x)| = in (|1 . 3j) and (|1.4p . This of course reflects the fact 
that inf spect(_ffs) = and therefore no exponential decay of the heat kernel is possible. 

On the other hand, Laptev and Weidl showed in |LWj that under certain conditions on B the 
operator Hb satisfies a Hardy type inequality 

Hb > r^ (1.5) 

1 + \x\ 2 

in the sense of quadratic forms on H 1 (M. 2 ), see also [WJ . Inequality (|1.5p implies that Hb is a 
subcritical operator. The criticality theory then suggests that the integral 

c 

e- tHB (x,y)dt (1.6) 

'o 

should be finite for all x ^ y. Hence in the limit of large times the magnetic heat kernel e~ tHB (x, y) 

should behave differently than the heat kernel of the usual Laplace operator in R 2 . Our motivation 

is to find out how exactly the large time behavior of e~ tHs (x, y) depends on the magnetic field. 

One of our main results, Theorem 14. 11 shows that for radially symmetric and weak magnetic fields 

the time decay of e~ tHn (x, y) is completely determined by the total flux of the magnetic field. The 

key point of the proof is to show that e~ tHn (x, y) is asymptotically (as t — > oo) equivalent to the heat 

kernel of certain two-dimensional Schrodinger operator with positive potential, see (|3.1|) . In section 

Owe establish some pointwise and L p — estimates on the magnetic semigroup e~ tHB in terms of the 

distance between the total flux and the set of integers, see Theorem 15.11 and Proposition 15.21 One 

of the main technical tools used in the proofs is Lemma 13.31 in which we derive a formula for the 

heat semigroup of certain family of one-dimensional Schrodinger operators. This also enables us to 

write an explicit expression for the heat kernel of the Aharonov-Bohm operator, Proposition [6TT1 and 

consequently to calculate its exact large time asymptotic, see Theorem 16.31 

Remark 1.1. Let us make a brief remark on the properties of the heat kernel under gauge transfor- 
mations. It is a matter of fact that the vector potential A is not uniquely determined by the mag- 
netic field B. However, if rotA = iotA = B G C(R 2 ,R), then there exists a scalar field such that 
A = A + V0. So the respective Hamiltonians Hb and Hb are unitarily equivalent; Hb — e 1 ^ Hb er 1 ^ , 
and their heat kernels are linked through the equation 



c 



tHB (x, y) = e 'W(*)-<^)) e - tHs (x, y). 



Hence changing the gauge does not change the time dependence of the heat kernel. In other words, 
the decay rate in time is gauge invariant. 

2. Preliminaries 

Given two functions f\, ji on a set f2 we will use the notation f\ ~ fa to indicate that there exist 
positive constants c, C such that the inequalities c fi < f% < C f\ hold on Vt. Accordingly, the 
notation fi(t,x) ~ f2(t,x) as t — > oo means that f\ ~ fi holds for all t large enough. Moreover, 
given two points x,y £ M 2 , we will often use the polar coordinate representation e~* B (x,y) = 
e~ tHB (r,r' ,6,6') of the heat kernel which corresponds to the identification x = r(cos8,smd) and 
y = r'(cos6' ,sm9'). Finally, we denote R.+ = (0,oo) and M + = [0, oo). We will need the following 
hypotheses. 

Assumption 2.1. Let the magnetic field be given as B(|x|), x € R 2 , where B : K + — > R is a 
continuous function with the support contained in some interval [0, R], R < oo. 

We introduce a vector potential A : R + x [0,27r) — >• R 2 which in polar coordinates (r,0) reads as 

follows 

1 f 
A{r,6) =o(r)(- sin0, cos 0), a{r) = - B(t)tdt. 

r Jo 
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Then A generates the magnetic field B. Hamiltonian Hb is associated with the closed quadratic form 

Qb[u}= / (\d r u\ 2 +r- 2 \id e u + b{r)u\ 2 )rdrd6, u £ iJ 1 (R+ x (0, 2tt)), (2.1) 

Jo Jo 

where 

b(r) = ra(r) = [ B{t)tdt = — / B(|x|)dx 

JO 2?r •/•(>|<r} 

is the flux of the magnetic field through the disc of radius r centered in the origin. Moreover, we 
denote by a the total flux of the magnetic field through the plane. By assumption 12.11 we have 

b(r)=a Vr>R. (2.2) 

By expanding a given function u £ L 2 (R + x (0, 2w)) into a Fourier series with respect to the basis 
{e"" e } me z of L 2 ((0,27r)), we obtain a direct sum decomposition 

L 2 (R 2 ) = X> £m ' (2.3) 

where £ m = {.g € L 2 (R 2 ) : g(x) = / (r) e imB a.e., J °° \f(r)\ 2 rdr < oo}. Since the magnetic field B 
is radial, the operator Hb can be decomposed accordingly to the direct sum 

H b= 51©(/* m ®id)n m , (2.4) 

meZ 

where /i m are operators generated by the closures, in L (R+, rdr), of the quadratic forms 
defined initially on Cq°(0, oo), and n m : L 2 (R 2 ) — > C m is the projector acting as 

M(M) = 7^ / " ^ m{6 - 6 ' ) u{r,e')de'. 

27r Jo 
Note that n m commutes with /i m (8> id. Hence the integral kernel of e~ fHB splits as follows: 

e- tHB (x,y) = i J] p m (r,r',t) e ^- e '). (2.6) 

Here p m (r,r',t) denotes e~ thm (r,r') which is real and positive for all m £ Z. The idea behind the 
proof of Theorem 14. II is to show that if the magnetic flux is small enough, then the large time behavior 
of e~ tHe (x, y) is determined by the contribution from m = in ()2.6[) . 



3. Heat kernel of the reduced operators 

In this section we will study the heat kernels p m (r, r' , t) . First we prove a result which allows us 
to quantify the large time behavior of Po(r, r' , t). To this end we consider an auxiliary Schrodinger 
operator 

A=-A + ^$- in L 2 (R 2 ). (3.1) 

The operator A can be defined in the usual way through the corresponding closed quadratic form 

QA[f\ = f (|V/| 2 + f ^P) dx, I £ D(Q A ) = H\«. 2 ). 
Jr2 V \x\ 2 J 
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Lemma 3.1. Assume that b(-) is not identically zero. Then there exists a positive radial function 
h G C 2 (M 2 ) such that Ah = 0. Moreover, any such function satisfies 

( IdH if n + 

Finally, there exist positive constants C and c such that the heat kernel of A admits for all x,y Gl 
and all t > the following estimate, 

e-^(x,y) ~ C ^y) e - c ^_ (33) 

Here we use, with a slight abuse of notation, the same symbol for the function h on M 2 and for its 
natural identification on K + . 

Proof. Since b(\x\) = a for \x\ > R, the spectrum of A coincides with the positive half-line [0, oo). 
Hence by the Allegretto-Piepcnbrink theorem, see e.g. |MPj . there exists a positive solution u to the 
equation Au = 0. Since the potential term 6 2 (|a;|)/|a;| 2 in A is Holder continuous, see assumption 
12.11 the elliptic regularity ensures that u E C 2 (K. 2 ). The radial function h given by 

h(\x\) = / u(\x\,0)d0, 
Jo 

then also satisfies Ah = and for \x\ > R we have 

h(x) = h(\x\)^a\x\^+b\x\-^, a^{) (3.4) 

h(x) = h{\x\)=c + d\\og\x\\, a = Q. (3.5) 

The positivity of h implies that a > 0, d > 0. On the other hand, h satisfies r(rh'(r))' — h(r)b 2 (r) 
with r = \x\ and therefore it is easy to see that h is an increasing function of r. This means that 
a > 0, d > 0. A straightforward verification now shows that the manifold K 2 equipped with the 
Lebesgue measure and the function h satisfy hypothesis of |Gr05| Thm.lO.lO.(i)]. The latter yields 
the heat kernel estimate (13.31). □ 



Corollary 3.2. There exists a positive radial function h G C 2 (IR 2 ) such that Hb h = 0. Moreover, 
if &(•) is not identically zero, then any such function satisfies (|3.2p . 



Proof. This follows from Lemma 13.11 and the fact that the operators Hb and A coincide on the set 
of radial functions. □ 

In order to control the terms in (12.6[) with m ^ we will make use of Lemma 13.31 below which gives 
an explicit formula for the heat semigroup generated by the operators which are associated with the 
quadratic form 

Qp[u]= f (u' 2 + ^ju 2 jrdr, [3 E R (3.6) 

defined on Cq°(M.+). This form is closable, see e.g. |Da| Sec. 1.8], and its closure generates a self- 
adjoint operator Tip in L 2 (R+, rdr). By the Beurling-Deny criteria Tip generates on L? (R+ , rdr) a 
symmetric submarkovian semigroup e^^* 3 . Let e~* w,3 (r, r') be its integral kernel. 

Lemma 3.3. Let Tip be the operator in L 2 (R+,r<ir) associated with closure of the form Qp. Then 
for all r, r' e M+ and all t > it holds 

1 (rr> 



e- tn ^r,r') = -l m [ — )e- r -^^, (3.7) 



where Lp\ is the modified Bessel function of the first kind, see e.g. [AS[ Chap. 9]. 
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Proof. Consider the operators 

L l3 = UH f sU- 1 in L 2 (R + ,o>), (3.8) 

where U : L 2 (R + ,rdr) — > L 2 (M. + ,dr) is a unitary mapping acting as (Uf)(r) — r 1 / 2 /(r). Note that 
Lp is subject to Dirichlet boundary condition at and that it coincides with the Friedrichs extension 
of the differential operator 

_l_ + l_A 

dr 2 r 2 

defined on Cq°(R+). Denote by D(Lp) the domain of Lp. Now let A be a complex number from 
some fixed neighborhood of R+ . A straightforward calculation using the standard technique of the 
Sturm-Liouville theory shows that the integral kernel of the resolvent operator (Lp — A)" 1 for r < r' 
is given as follows 

(Lp - X)- l (r,r') = y\^ J\p\(rVX) (j m (r'VX) +iY m (r'V\)) , ImA > 

(Lp - X)-\r,r') = -^VrV J\p\(rVX) ( J\p\(r'VX) - iY m (r'V\)) , ImA < 0, 

where J\p\ and Y\p\ are the Bessel functions of the first and second kind respectively. Next wc 
introduce the function g(r, A) = y/r J\p\(r\/X) , and note that Lp g — Xg and g(0, A) = 0. Hence the 
Weyl-Titchmarsh-Kodaira Theorem, see [DSchl Chap. 13], says that 

W fi LpWp 1 ip(X) = Xtp(X) ) v^Wp{D(Lp)), (3.9) 



where the mapping Wp and its inverse Wo 1 given by 



(Wpu){X)= u(r)Vr~J m (rVX)dr, (WCV)W = / <pWV* J\p\(rVX) — (3.10) 



dX 



ii 



2 



defined initially on Cg°(R+) extend to unitary operators from L 2 (R + ) onto itself. Given / £ Cq°(R+), 
in view of (|3.9|) we then get 

(e-*^ /) (r) = (W^ 1 e~ tx Wp /) (r) = J VrV J e~ tx J w (rVX) J m (r'Vx) — f(r') dr' 



1 r°° — / rr '\ r '+r» 

2t ./„ ^^'U 



rr' I\p\[ — )e- r -^ f(r')dr', (3.11) 



where we have used Fubini's theorem to switch the order of integration and [Erdl Eq.8.11(23)] to 
evaluate the A— integral. Moreover, since Vr 7 Lpi(rr' /2t) e « £ L 2 (R+) for all r, t > 0, see |AS| 
Chap. 9. 7], identity (J3.11JI extends by density to all / £ L 2 (R + ). Hence 

e- tLl '{r,r r )'j=—y/rr'I\B\[ — )e — ~ (3.12) 



2t " l/31 V 2t 
is the integral kernel of e~ tLf> , and by (|3.8|) we conclude that 



-™?(r,r') = -=L e- tL »(rS) = i I m (jp) e" 1 ^. (3.13) 



D 



Lemma 3.4. Let \a\ < 1. Tften /or a?/ ijel 2 it holds 



lim e-^^^bod^Myl,*))-^!. (3.14) 
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Proof. Operator A admits the decomposition 

A=^2®(A m ®id)n m , (3.15) 



1116Z 
1( 



where A m are operators in L (K + , r dr) generated by the closures of the quadratic forms 



Jo \ ' / 

defined on Cq°(0,oo). Note that Aq = ho and hence 

e- fA (x,y)=p (\x\,\y\,t)+Yl e~ tA -{\x\,\y\)e im ^- e '\ (3.16) 

In order to estimate the sum on the right hand side of the last equation, we note that by the Trotter 
product formula 

e~ tAm {r,r') < e~ mm (r,r') Vr,r'eK+, Vm^O, (3.17) 

where T-L m is the operator defined in Lemma 13.31 By the same Lemma we get 



| Y, e-^(M,b|)e^-*')| < Cz ]T V|(z), z:=^, (3.18) 

where the constant C depends on x and y. Assume first that a^O. From the integral representation 

Iu(z)= ; Z \, 7TT / (1 - s 2 ) v -? e zs ds (3.19) 

vy ' 2T(i/ + |)r(|) J-i ' y ' 

for J„, see e.g. |AS( Chap. 9], it is then easy to see that 

_n — |q 

limsup i 1+ l"l | ]T e- Mm (kl, |y|) e m{e ~ e ) | < c limsu P E on TV 7^ 

^°° m^O 2 ^° n>l ^ U?1+ 2 J 

< c limsup z 1Hq| V ^ r- = 0, 

n>l v *■' 

where c' depends on x and y. Since 

e~ M (x ;2 /) ~ r 1 '^ t^oo, a^O, 

by Lemma [3.11 we conclude from (|3.16p that equation Q3.14p holds true in the case a ^ 0. On the 
other hand, if a = 0, then Lemma 13 . 1 1 gives 

e~ tA {x,y) ~ i _1 (logf) -2 t -> oo, a = 0. 

From (J3TT5|) and (|3"7P9"f we find 

(logz) 2 2;™ 



limsup t (log t) 2 | J^ e _ * A "(|a;|,|v|)e 4m(e - fl ' ) | < c limsup ^ 



.-n rr[ 2"F(n+|) 



This proves (pTi]) for a = 0. □ 

4. Large time asymptotic of e~ tHB (x, y) 

Below we formulate our main result regarding the large time behavior of the magnetic heat kernel 
e -tH B ^ yy n snows that if the magnetic field is sufficiently small, then the decay rate of e~ tHB (x, y) 
is completely determined by the total flux a. 
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Theorem 4.1. Let B(x) satisfy assumption 1 2. 1\ and suppose that \b{r)\ < 1/2 for all r £ R+. Let 
h G C 2 (R 2 ) be a positive radial function such that Hsh = 0. Then there exist constants C and c 
such that the inequalities 

c < liminf t 1+ l«l e " tg " ( ^f ) < limsup £ 1+H ^"^'^ < C, a^O (4.1) 

t-voo n{x)h{y) t^oo h{x)h{y) 

and 

c < lim inf t (log tf e ~'" B ^ ] < lim sup t (log i) 2 ^w'^ < C, a = (4.2) 

t^oo n{x)h{y) t-+oo h{x)h{y) 

hold true for all x,y £l 2 . 

Remark 4.2. Similar connection between the large time asymptotic of the heat kernel e~ tp (x,y) 
and the ground state of the corresponding generator is known when P has an eigenvalue at the 
bottom of its spectrum, see e.g. |CK1 [P |ISi3j . 

Remark 4.3. Equation (|4.2|) shows that e~ tHn (x,y) is integrable with respect to t at infinity even 
if the total flux is zero. The latter reflects the fact that Hb satisfies a Hardy type inequality also in 
this case, see |W] . 

Proof of Theorem \4-l\ The existence of the ground state h is guaranteed by Corollarv l3.2l By Lemma 
13.11 it suffices to show that 

lim e J H ^ ,V \ =1 Vx,2/eR 2 . (4.3) 

t-too e ttiB (x, y) 

Let a ^ 0. By assumption we have (b(r) + m) 2 > to 2 /4 for all m ^ and all r G R+. Hence the 
Trotter product formula gives 

p m (r,r',t) < e _tWm / 2 (r,r') Vr,r'eR+, t > 0, m ^ 0. 

With the notation of equation (|3.18[) we get from (J3.19I) and Lemma 13.31 

limsup i 1+ H J2 Pm(r,r',t) < c limsup z £ I |m/2| (z) < c lim £ ^ ^ - = 0, 

m^O m^O n>l 

where we have used the fact that \a\ < 1/2. In view of equations (|2.6[) . (|3.3|) and Lemma f3.4[ this 
proves (|4.3p . If a = 0, we obtain in the same way as above 

limsup t (log t) 2 ^2 Pm(r,r',t) = 0. 

t— s-oo , 

Equation (|4.3p thus holds also in this case. □ 

In the case \a\ > 1/2 we give an asymptotic upper bound on the heat kernel. 

Proposition 4.4. Let B{x) satisfy assumvtion \2. 1[ Let g — minfcgz \k + a\ be the distance between 
the flux a and the set of integers. Then there exists a constant C such that 

limsupt 1+e |e- tHB (x,2/)| < C(l + \x\) e (1 + \y\) e (4.4) 

t— s-oo 

holds for all x,y £ R 2 . 

Proof. We introduce the operators T m generated by the quadratic forms 

defined initially on Co(R+) and then closed in L 2 (M.r dr). Here 0(-) denotes the Heaviside function. 
By the Trotter product formula we have 

e - th ™(r,r') < e- tTm (r,r') Vr,r'eR+, m e Z. (4.5) 
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In view of (|2.2j) it follows that the functions ip m £ C 1 (R+,R+), defined by 

1 /r\ CT ™ 1 / R\ am 
i> m {r) = 1, r<R, VmW = - [-) + g I 7 ) r - R > ° m = \ a + m \ ( 4 - 6 ) 

solve the Cauchy problems 

Q(r _ R) (Hr)+m) 2 ^ = (r ^y ) ^ m{R) = ^ ^ (i?) = Q 

The operators 

5*™ = ipm T m ip m in L 2 (M+, ipm(r) rdr), 
are thus unitarily equivalent to T m and their heat kernels satisfy 

e -«- m(r y) =^(r)^(r') e - iS "( r ,r'). (4.7) 

A direct calculation shows that 5Vn is associated with the quadratic form 

/>00 

s m [u] = t m [u VVn] = / (u') 2 ip^r dr, u £ D(s m ) = H 1 (R+ , ^ r dr) . 



We now apply Theorem I A. II with /i(x) = z/(x) = :c0 m (a:), p = 2 and g = (2 + 2a m )/a m . Hence for 
each m there exists a constant c m , such that 



s m [u]>c m { u °™ tp m rdr) Vu£D(s m ). (4.8) 

By the Beurling-Deny criteria, S m generates on L 2 (WL+,ip m (r)dr) a symmetric submarkovian semi- 
group e~ tSm . This allows us to apply [Dai Thm.2.4.2], see also [Var j . to obtain 

||e^ s "1oo, 2 < C m t- 1 - ± ^ 

for some constant C m . By duality this implies that 

sup e~ ts ™(r,r>) = ||e- tS '"||oo,i < IK' 5 '"!!^ < C?*- 1 -*" V* > 0. 

rr i 

In view of equations (J4.5I) , (I4.6P and (I4.7[) this yields 

e"* Tm (r,r') < C^ t- 1 - CT -(l+r) CT '"(l + r')' Jm i > 0, r,r'eM+. (4.9) 

Now define uq :— mi{n £ N : n > 2sup r>0 |fr( r )l}- From (I4.5[) and (|4.9[) we obtain 

limsupi 1+e V p m {r,r',t) < C(l + r) e (l + r') e . (4.10) 

m— -no 

To estimate the rest of the sum in (12.6[) we note that 

o 
772. 

(6(r) + m) > V r > 0, m : |m| > no- 

Hence mimicking the arguments used in the proof of Theorem 14.11 it is easy to see that 

limsupi 1+<? V p m (r,r',t)=0. 

\m\>n 

By (|2.6|) this completes the proof. □ 
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5. Heat kernel estimates 



In this section we use Theorem 14.11 and Proposition 14.41 in order to prove certain point-wise heat 
kernel estimates. We use the notation introduced in Proposition 14.41 i.e. g — min/jgz \k + a\. 



Theorem 5.1. Under assumvtion \2.1\ there exists a constant C such that the inequality 

e- tHB (x, x) < C min {t~\ (1 + \x\) 2e i" 1 " 5 } (5.1) 

holds for all igt 2 and all t > 0. 
Proof. Adopting the notation of the proof of Theorem 14.41 it follows from (|4.9[) that 

Pm(r, r, t) < CT 1 -^ (1 + r) 2 ^ Vm e {-no, ■■-, no}. 
On the other hand, the diamagnetic inequality 

\e- tHB (x,y)\ < e tA (x,y) = ^- t e- 1 ^ 1 , x,yeR 2 , t > 0, (5.2) 

see e.g. |AHS[lHSllSi2| . clearly implies that p m (r,r,t) < l/(2t) for all m. Hence 

p m (r,r,t) <C(r 1 -« r ~(H-r) 2 ' r »)^rt-f 1 -=fc- ) = C T 1 ' 8 {I + rf e \m\ < n . (5.3) 

I i2 

Next wc introduce the variable z — ^-. From the proof of Theorem 14. II arid from Lemma 13.31 we get 



J2 tp m (\x\,\x\,t)< z-« £ te- tn «"*(\x\,\x\) = cz-°e-' £ I\ m/2 \(z). (5.4) 

|m|>no |m|>no |m|>no 



On the other hand, inequality (|5.2|) shows that 

J2 tp m (\x\,\x\,t) < Y,tPm(\x\,\x\,t) =2irte- tHB {x,x) < \. 



|m|>no m£2 

This in combination with (|5.4|) gives 



sup—— V p m (\x\,\x\,t) < cmax sup z~ e V I\ m/2 \(z) , sup z~ e \ < C. 

t,r>0 t y , i L x<l , i z>l J 

\m\>iiQ — |m|>no 

Indeed, in view of p. 191) the series J2\ m \>n I\m/2\( z ) converges uniformly with respect to z on [0, 1]. 
Hence ;>H Q J2\ m \>n I\m/2\i z ) is continuous on (0,1] and since it tends to zero as z — > 0, see the 
proof of Theorem 14.11 it is bounded. From equation (|5.3p we thus get 

Y,Pra{\x\,\x\,t) <C{l + \x\f e t' 1 - 8 , Viel 2 , £>0. 

mGZ 

The statement now follows by (|2.6[) and (|5.2[) . D 



As a consequence of inequality (|5.1I) we get an estimate on the norm of e B acting on certain 
weighted LP spaces. To formulate our result we introduce the following family of subspaces: 

\_ 

L$ := {/ : ||/IU ? < ™} , 11/11*5 := (^ \f\ P (1 + \x\f <*c) * , [3 G R. 

We then have 

Proposition 5.2. Let as sumptions ] 2. 1\ be satisfied. Assume that p € [1,2] and Zei g € [2, oo] be such 
that - + - = 1. TTien for any (3 > 2 + 2 g there exists a constant C = C(g, (3) such that 

\\e- tHB \\ L$ -, L , < Cr?~i Vi>l. (5.5) 



\ e -tH B 
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Proof. We use the shorthand kit, x, y) = e~ fHB (x, y). Since e~ tHB is self-adjoint, we have k(t, y, x) — 
k(t,x,y). The semigroup property of e~ tHs and the Cauchy-Schwarz inequality then yield 

\k(2t,x,y)\= / k(t,x,z) k(t,z,y)dz <( / \k(t, x, z)\ 2 dz ) I / \k(t, z, y)\ dz) 

Jr 2 \Jr 2 ' v Jr 2 ' 

= y/k(2t,x,x) y/k{2t,y,y). (5.6) 

This in combination with estimate (|5.ip and diamagnetic inequality (|5.2p gives 

\k(t, x,y)\< CT 1 -?* 1 ^ (1 + |x|r e (l + |y|r e , V M i, M2 G [0, 1]. (5.7) 

Now fix / e L 2 (R 2 ) and let t > 1. Chose m = /j, 2 = 1 in ([5J]) . In view of (|5T6| . Cauchy-Schwarz 
inequality and Fubuni's theorem we have 

f\\h= f I k(t,x,y)f(y)dy\ 2 dx < \\f\\% f f \k(t,X,y)\ 2 (l + \y\)^ dydx 

JR 2 JR 2 ? JR 2 JR 2 

= ll/ll 2 / k(2t,v,y)(l + \v\)-f>dy < Cf 1 "" ||/||i ? . (5.8) 

This shows that 

\\e- tHB \\ L 2 ^L 2 < C'^ (1+e)/2 . (5.9) 

On the other hand, choosing fix = and /ia = 1 m (|5.7p it is easily seen that 

l|e-* HB |U^~ < Ct- 1 -^ 2 . (5.10) 

Inequality (|5.5[) now follows from (J5.9I) . (|5.10|) and the Riesz-Thorin interpolation theorem. D 

Remark 5.3. In the absence of magnetic field we have 

l|e~ tHo IU^ = l|e* A || i? ->£. ^CH 1 V/3 > 2. (5.11) 

Indeed, the upper bound in (|5. 11|) follows by mimicking the proof of Proposition 15.21 with k(t,x,y) 

, ,2 

replaced by e tA (x,y) = e J* /(47rt). This leads to equations (|5.9p and (|5.10[) with g = 0. In order 
to prove the lower bound in (|5.11[) let us consider the solution of the heat equation with the initial 
data f(x) = e~' x '". An easy calculation gives 

u(t,x) = (e tA f) (x) = -^ e-±&, ||u(V)IU« = 0(1 + 4*)*?*. 

Proposition 15.21 thus says that the L q norm of the solution to the heat equation 

d t u + Hbu = 0, w(0, a;) = f(x), 

decays faster (with respect to the case B = 0), if we restrict the initial data / to a smaller subspace 
of L P (R 2 ). Note also that similar estimates were recently obtained, in the case p = q = 2, for the 
heat semigroup of Dirichlet-Laplace operator in twisted waveguides; see |KZj . 

6. Example: The Aharonov-Bohm operator 
A natural question which arises from theorem 14. II is whether the limit 

lim t 1+M e- tHB {x,y) (6.1) 

t— >oo 

always exists and how it depends on x and y. In this section we calculate the limit (j6.ll) in the case 
of the so-called Aharonov-Bohm magnetic field. This field is characterized by the property that the 
flux b[r) through a disc of radius r is constant. It is generated by the vector potential A whose radial 
and azimuthal components (in the polar coordinates) are given by 

A(r,0) = (ai(r,e),a 2 (r)), a x = 0, a 2 (r) = (o,"). (6.2) 



LARGE TIME BEHAVIOR OF MAGNETIC HEAT KERNELS 11 

The associated operator (iV + A) 2 defined on C£°(R 2 \{0}) has deficiency indices (2, 2), see [ATJlPR]. 
We will consider the Hamiltonian H a as its Friedrichs extension. In other words, we define H a as a 
non negative self-adjoint operator in L 2 (R 2 ) generated by the closure of the quadratic form 

/■2-7T /"OO 

QM = / (\d r u\ 2 + r- 2 \(-id e + a)u\ 2 )rdrd9, u G C o °°((0, oo) x [0, 2tt)). 

Jo Jo 

Proposition 6.1. Lett > andr,r' £ R + . Then the heat kernel of the Aharonov-B ohm Hamiltonian 
H a is given by the absolutely convergent series 

Proof. We note that 

H * = J2 ® ( w ™+« ® id ) n '»' ( 6 - 4 ) 

where H m + a are the operators in L 2 (R + ,rdr) defined in Lemma \3. 31 Hence 



(*»») 4^ e -^+=( r y )e -(^') , (6.5) 



e -tK, 

' 2vr 



Equation (|6.3|) now follows from Lemma 13.31 The absolute convergence of the series is easily seen 
from the integral representation of the Bessel function /„, see equation (|3.19|) . □ 

Remark 6.2. For a £ Z we get by (Ml Eq.9.6.33] 



e 



* H -( x ,y) = — e-^e^^-^e^ '-^ = — e"^ e ^' -*» . (6.6) 

V ' y; Ant Ant V ; 

This reflects the well known fact that for integer values of the flux the Aharonov-Bohm operator 

is unitarily equivalent to the Laplacian in L 2 (R 2 ) under the unitary mapping / i-> e~ ta6 f , see 

also Remark 11.11 Equation (|6.3|) also implies that it is no loss of generality if we suppose that 

a £ [-1/2, 1/2]. 



1 / rr '\H 



Theorem 6.3. We have 

lim ^H e-» (*, y) = ^^^ (^ j */ |a| < 1/2, (6.7) 

1 

t ^e-^( X ,y) = 1 J m {^-y ( l + e -(-')) if a = ±1/2. (6.8) 

Proof. From equation (|3.7[) and the asymptotic behavior of l v (z) for small z, see [AS! Chap. 9], we 
get 

\m+a\ r i\m+a\ 

lim t 1+lm+al e -**™+« (r, r') = „, ' , / : . (6.9) 

*^oo v ' ; 2 2 l™+ Q l+ 1 r(l + |m + a|) V ; 

Assume first that \a\ < 1/2. In view of (|3.19p we obtain 



i rr'\ ( rr ')\ m + a \ _, r _ r ' ) 2 

2t J - 2l m + Q lr(|m + a| + i)r(i) 



Since inf m ^o \fn + cx\ > \a\, it follows that 

Hm *W e -^ J] 7|„ l+Q | (^) e ""«^'> =0, 






which, in combination with (|6.9p . proves equation (|6.7|) . The proof in the case \a\ = 1/2 follows the 
same line. □ 
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Appendix A. 



For the reader's convenience, and also because equation (13) of [Ml Sec. 1.3.1] contains a missprint, 
we recall below a simplified version of [M, Thm.1.3.1.3]. 

Theorem A.l (Maz'ya). Let 1 < p < q < oo and let /i, v £ L 1 (R + ) be nonnegative. Then the 
inequality 

(f \u\ q fi(x) dxY < C ( f \u'\ p p(x)dxY (A.l) 

holds for all u £ W 1 ' P (J& + , v(x)dx) and some constant C , independent of u, if and only if 

sup I / fj,(x)dx\ I / u(x)^'p TrT dx I < oo. (A. 2) 
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